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I. INTRODUCTION 

Usually, classical spin models have global symmetries, 
that is, they are invariant under a transformation which 
is the same at all points in spacetime. Instances of such 



theories are Newtonian dynamics under Galilean trans- 
formations, special relativity under Lorentz transforma- 
tions, or the choice of the zero point of energies in any 
system. These symmetries are motivated by the homo- 
geneity and isotropy of space, and the homogeneity of 
time, which render the choice of the origin of the co- 
ordinate system, as well as the orientation of its axis, 
arbitrary. The same sort of symmetries apply for usual 
statistical models, such as the Ising or Potts model, where 
a global rotation of the classical spins leaves the physics 
of the system invariant PQ. We will henceforth refer to 
these models as Standard Statistical Models (SSMs). 

The notion of global symmetry can be lifted to the 
notion of local symmetry, where the applied transfor- 
mation is point-dependent. This is precisely the case 
for gauge theories, which describe the most fundamen- 
tal interactions in nature, like quantum electrodynamics 
(QED), weak interactions, and quantum chromodynam- 
ics (QCD). These symmetries are physically motivated 
by the conservation of certain quantities at every point 
in spacetime. More precisely, the conservation of electric 
charge, the weak neutral current, and the color charge 
induce the U(l), SU(2), and SU(3) gauge symmetries in 
QED, weak interaction theory, and QCD, respectively. 

Lattice gauge theories (LGTs) are lattice formulations 
of gauge theories [2]. That is, they are theories in which 
the variables sit at the edges of a lattice, and the Hamil- 
ton function of the system exhibits some local symme- 
tries. As we will elaborate below, they are interesting 
not only as a means to study gauge theories, but also as 
a new class of statistical models per se. 

Given this variety of models (SSMs, with global sym- 
metries, and LGTs with local symmetries), one may won- 
der if it is possible to relate them, that is, to find a map- 
ping such that having information of one model automat- 
ically yields insight into another model. The purpose of 
this paper is to show that a general mapping of this kind 
exists. More precisely, we will prove that the partition 
function of any (Abelian, discrete) classical spin model 
can be expressed as the partition function of one specific 
model, which is enlarged and which contains inhomoge- 
neous coupling strengths. Note that the partition func- 
tion is the crucial quantity of a model, since, for SSMs, 
one can obtain all thermodynamical properties by taking 
derivatives thereof [3], and for LGTs, one can compute 
relevant quantities such as Wilson loops [IHS] . Following 
Ref. [7J, we will call this a completeness result, since it 
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shows that the partition function of one specific model 
can specialize (by being defined on a large enough lattice, 
and by tuning in it the appropriate coupling strengths) 
to the partition function of any other (Abelian, discrete) 
classical spin model. 

Similar completeness results have been proven before. 
A very general result was proven in Ref. [7], where it 
was shown that the 2D Ising model with magnetic fields 
is complete for all other (Abelian, discrete) classical spin 
models. However, the complete model needed to allow for 
complex coupling strengths in order to specialize to other 
models, and thereby hindered a physical interpretation 
of the complete model. In Ref. [S] several completeness 
results with real parameters were obtained, and it was 
shown, e.g. that the 3D Ising model is complete with real 
parameters for all other Ising models with fields. While 
similar results for other types of models were shown, none 
was as general as the original result of [7J , thereby giving 
the impression that there was a tradeoff between the real 
parameters and the generality of the completeness re- 
sults. This tradeoff has proven to be false, for we showed 
in Ref. [9 that one can obtain a complete model for all 
others (Abelian, discrete) classical spin models employ- 
ing only real parameters. Moreover, we showed that the 
Hamilton function of a subsystem of the complete model 
equals the Hamilton function of any given classical spin 
model. 

One may wonder why LGTs are interesting from the 
point of view of the completeness results. The reason 
is that gauge theories are models with local symmetries, 
which are a generalization of systems with global symme- 
tries. Thus, there are chances to obtain more powerful 
completeness results than with SSMs, and this is pre- 
cisely what we have found. 

In this work, we elaborate on Ref. [5] and illustrate that 
result with specific examples and constructions. More 
precisely, in [9] it was proven that the partition func- 
tion of any Abelian, discrete classical spin model can be a 
mapped onto the partition function of a lattice gauge the- 
ory on a 4~ dimensional square lattice with gauge group 
Z 2 (the 4D Z 2 LGT). Here, we go beyond Ref. [5] in the 
following points: 

• We present a quantum formulation of the partition 
function, where the latter is expressed as the scalar 



product between two states (Sec. Ill A I. Here, we 
show that one of these two states is a stabilizer 
state, which reveals some symmetries of the parti- 
tion function; 

• We present a detailed analysis of the merge and 
deletion rule (Sec. 



IIIB), the basic tools used to 



manipulate and transform one interaction pattern 
to another; 

• We present a rigorous proof of the method to obtain 
general n— body interactions (Sec. IIIC); 



teractions are present (termed fc-cliques), for all 
k = 1, . . . , n; in particular, we show how particles 
propagate and their paths turn around in the lat- 
tice, how one can construct 4- and 5-body Ising- 
type interactions, and the explicit layout of all k- 
cliques in the 4D lattice (Sec. HID); 



• In Sec. HIE we give a detailed explanation of the 
class of models embraced by the main result of this 
work; 

• We provide a new discussion of the main result in 
terms of the space of all theories (Sec. Ill H I . 



Moreover, we will explore some of the consequences of 
this result, going beyond [H| in the following aspects: 

• In Sec. |IV| we give a general observation concerning 
all the completeness results; 

• We illustrate our result by explicitly computing the 
Wilson loop of a lattice gauge theory, and the mag- 
netization of the 2D Ising model from the partition 
function of the 4D Z 2 LGT (Sec. 0. We also dis- 
cuss how one can map models with global symme- 
tries to a model with local symmetries. We also give 
an optimized construction of the 2D Ising model. 



• In Sec. VI B we provide an explicit construction of 
the 4-clique required to compute the mean-field 
theory. 

This paper is structured as follows. First, we will 
present a brief introduction on lattice gauge theories for 
the non-expert reader (Sec.|ll]). In Sec. Ill we will prove 
our main result, namely that the partition function of any 
Abelian discrete SSM and any Abelian discrete LGT can 
be mapped to the partition function of the 4D Z 2 LGT. 
In Sec. |IV| we will make a general observation on how 
to compute observables of a model from the partition 
function of the complete model. Then we will explore 
some applications of this result in Sec. [V] In Sec. |VI| we 
will focus on two implications of the main result, namely 
the computational complexity of the 4D Z 2 LGT and a 



We provide the explicit construction of interaction 
patterns where all possible /c-body Ising-type in- 



new mean-field theory for Z 2 LGTs. Finally, in Sec. VII 
we will present some extensions of the main result, that 
is, we will show the completeness of other models. In 
Sec. |VIlT| we will summarize our results, discuss a possi- 
ble reduction of our result to 3D and give an outlook on 
the subject. 



II. BACKGROUND ON LATTICE GAUGE 
THEORIES 

LGTs were originally introduced by Wegner with the 
goal of obtaining statistical models with a vanishing mag- 
netization in all phases, but nonetheless with a non trivial 
phase diagram [10] , In this approach, LGTs are seen as 
statistical models with local symmetries, which exhibit 
novel features compared to SSMs. For example, because 
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local symmetries cannot be broken by local order pa- 
rameters (such as the magnetization) [11) . one needs to 
define non-local order parameters to witness the differ- 
ent phases of these models. An example of a global order 
parameter is the Wilson loop, which is a product of spins 
sitting at the edges (or, more generally, of gauge fields) 
over a closed path: W(C) — Yiiec Si - 

The Wilson loop was introduced to distinguish con- 
fined from deconfined (or Higgs) phases, since in the for- 
mer phase it exhibits an area law, whereas in the latter 
it exhibits a perimeter law |12| . It is the order parameter 
for LGTs, and it plays a similar role as the magneti- 
zation for SSMs. Another important quantity is the 't 
Hooft loop, which is related to the Wilson loop via a cer- 
tain duality transformation. Physically, the Wilson loop 
is an order variable, while the 't Hooft loop is a disorder 
variable, generalizing concepts similar to those in SSMs, 
which arise in the context of the Kramers- Wannier du- 
ality. In the Hamiltonian formulation of an LGT |13j . 
Wilson loops on a closed curve C, W(C), and 't Hooft 
loops on a closed curve C", W(C), are operators that sat- 
isfy a loop algebra which is a type of Weyl commutation 
relation, like for momentum and position operators: 

W(C)W(C) = z L{c > c ">W{C')W{C) (1) 

where L(C, C) is the linking number of the two loops 
and z is an element (phase) of the group Z g , which is the 
center of the gauge group. With this algebra it is possible 
to draw conclusions about the perimeter and area laws 
for both types of loops, and use them to characterize 
phases in LGTs. 

The notion of gapped versus gapless phases is also 
present in LGTs, but gapped phases have a richer struc- 
ture than in SSMs. In fact, both confined and Higgs 
phases are gapped, and this is why the Wilson loop is 
needed to distinguish them. 

A correlation length £, and thus a gap, can be defined 
by means of a two-point correlation function G c (/i,/2) 
between elementary faces fi and separated a distance 
I apart, as follows: 

G c (l):=(s h s h )-(s fl )(s h ) (2) 

where the subscript stands for the connected component 
of the Green function. Thus, when the system is gapped, 
we find a behaviour G c (l) ~ e~ 1 ^, typical of a system 
with a finite correlation length. The system is gapless 
when the correlation length becomes infinite. 

Now, with the aid of the correlation length and the 
Wilson and 't Hooft loops, we can give a more precise 
picture of the more relevant phases in an LGT [T5J [T5] 

Confined phase: it is gapped and the Wilson loop obeys 
an area law, while the 't Hooft loop follows a 
perimeter law. 

Deconfined (Higgs) phase: it is gapped and the Wilson 
loop obeys a perimeter law, while the 't Hooft loop 
follows an area law. 



Coulomb phase: it is gapless and both loops obey a 
perimeter law. 

In another approach to LGTs, Wilson independently 
introduced a more general class of LGTs, namely non- 
Abelian LGTs (as will be defined below), as gauge the- 
ories formulated on a discrete spacetime |12) . In this 
context, LGTs are cutoff regulations of gauge theories 
of strongly interacting particles. Gauge theories them- 
selves describe three of the four fundamental interactions 
in nature: electromagnetic, weak and strong, with gauge 
group U(l), SU(2) and SU(3), respectively. Only grav- 
ity has evaded a consistent quantum gauge formulation. 
Thus, LGTs offer a new way to tackle hard problems 
in gauge theories, either by analytical or by numerical 
calculations. The most prominent example is quark con- 
finement, which has been shown to exist in a 4D SU(3) 
LGT, but it is still unclear whether it survives in the 
continuum limit |16j . 

Generally, LGTs can be classified according to their 
gauge group: Abelian discrete LGTs, with gauge group 
Z 9 , Abelian continuous LGTs, with gauge group U(l), 
non-Abelian discrete LGTs, with the permutation group 
S„, and non-Abelian continuous LGTs, with gauge group 
SU(n). These theories have applications in a variety of 
fields of physics. We have already mentioned applications 
of continuous LGTs in the Standard Model, when one 
takes the limit of continuous spacetime. Non-Abelian 
discrete LGTs, with the permutation group S3, are used 
to describe magnetic monopoles and topological defects. 
In this work we will focus on Abelian discrete LGTs. 

Despite being the simplest instances of LGTs, Abelian 
discrete LGTs already have numerous applications in 
physics. Already the Z2 LGT on a 3-dimensional square 
lattice (3D Z 2 LGT) exhibits a nontrivial phase dia- 
gram with a confined and a deconfined phase (which is 
known from its duality relation to the 3D Ising model), 
thereby mimicking the confinement properties that ap- 
pear in strong interactions, described by SU(3) LGTs. 
On the other hand, Z g LGTs can also be used as toy 
models to approximate U(l) LGTs. To do so, one only 
needs to let the number of levels of every particle q tend 
to infinity, as we will do in this work. A third applica- 
tion was introduced by 't Hooft, who observed that Z q 
is the center of SU(n) [14] . Thus, gaining insight into 
Z q LGTs may shed light on SU(n) LGTs, and, as men- 
tioned above, the latter are used to describe weak and 
strong interactions. Still another application is the fact 
that the hamiltonian formulation of Z 9 LGTs in d di- 
mensions can be mapped to a quantum theory in d — 1 
dimensions [13 . Furthermore, Abelian discrete LGTs are 
studied in quantum error correction of topological quan- 
tum memories, since a faulty syndrome measurement of 
these memories can be mapped to a Z g LGT with ran- 
domness [T7] . Finally, Z g LGTs can also be used to study 
spin glasses [15] . 

Formal definition of an Abelian discrete LGT. We con- 
sider a standard definition of an Abelian discrete LGT 
using a Wilson Hamilton function in terms of face inter- 
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actions, with gauge group Z q . That is, classical spins 
sit at the edges e € E of a d— dimensional lattice, and 
they have q levels, s e = 0, 1, . . . , q — 1. They interact via 
the faces / € F of the lattice; thus, a face of A: sides re- 
flects a fc-body interaction. More precisely, the Hamilton 
function of the system reads 



H(s) 



E J / Re 



fEF 



eedf 



(3) 



where J/ is the interaction strength at face /, Re stands 
for the real part of the expression, and e € df refers to the 
spins at the boundary of face /. Here s := (si, . . . , si e i) 
stands for the spin configuration of the system. For 
q = 2 (and any k and d) we will refer to these inter- 
actions as "Ising-type interactions" , and to the model as 
d-dimensional Z 2 LGT (dD Z 2 LGT) (by default be- 
ing defined on a square lattice, i.e. k = 4). Notice 
that in this case each face interaction takes the form 
J\...k{— l) Sl+ ■■■ +Sk , that is, it only depends on the par- 
ity of the k adjacent spins. Note also that the Hamilton 
function ^ corresponds to a "pure" lattice gauge the- 
ory (the quotations are due to this concept being usually 
defined for U(l) LGTs or SU(n) LGTs), since we only 
have "gauge fields" at the edges and there are no particles 
("matter fields") at the vertices. 

The Hamilton function (|3| is invariant under the gauge 
transformation 



9v = J| X e 

e:eadj v 



(4) 



where eadj v are the edges adjacent to a vertex v, and X e 
is defined as X e : s e i-» (s e + l) m od q (i-e. it is the classical 
analogue of the generalized Pauli operator X) . The gauge 
group is generated by these transformations applied on 
any vertex v, i.e. Z g :— {g v ,^v G V). LGTs are defined 
on an oriented lattice, that is, a lattice in which each edge 
has a tail (one of its end vertices) and a head (the other 
end vertex) . Then one applies the gauge operation to an 
edge adjacent to v if v is the tail of that edge, and the 
inverse gauge operation if v is the head of edge, where 
this convention is arbitrary. In this manner, every closed 
face (i.e. where all edges are oriented either clockwise or 
counterclockwise) will be gauge invariant, since a gauge 
operation in any of its vertices will cancel. In fact, one 
can assign a different orientation to the same edge for 
each face where it participates (see, e.g., [5] p. 694), since 
it is only relevant that every face is closed. Note that only 
in 7Li LGTs it is irrelevant to have an oriented lattice, 
since every element in the gauge group corresponds to 
its own inverse. 

In our proof of the main result, we will make use of 
the gauge fixing of edges. This consists of fixing spin 
values to zero at the expense of reducing the gauge sym- 
metry of the model. The resulting model is physically 
equivalent to the original one as long as one fixes spins 
at edges forming at most a maximal tree (i.e. not form- 



ing a closed loop [H]). In other words, all these mod- 
els would be described by the same "effective" Hamilton 
function (i.e. acting on the actual degrees of freedom, 
which arc those described by the Hamilton function con- 
strained by its symmetries), and in this sense they can 
be regarded as belonging to the same equivalence class. 
However, if edges fixed by the gauge form a (closed) loop, 
the physics described by the model change. Intuitively, 
this can be understood in several ways. First, if we fix 
all spins around the face, then the interaction in this face 
is fixed (Jf(— 1)°), so this is equivalent to deleting that 
face. Hence, this would amount to creating a "hole" in 
that face or changing the topology of the lattice. An- 
other way to see it, is that in lattice gauge theories order 
parameters are products of gauge fields around a closed 
loop. Hence, one cannot fix all these gauge without af- 
fecting the order parameter. 

The partition function of an Abelian discrete LGT is 
defined as 



J LGT 



(5) 



where /? = l/(kgT) is the inverse temperature. Al- 
though the partition function has no physical interpreta- 
tion as such, it is important because one can compute all 
other thermodynamical quantities as a function of it [3] . 
For example, the mean energy is obtained by taking the 
derivative with respect to (3 of the logarithm of Z 



we will give examples of these calculations in Sec. IV 



Thus, complete knowledge of the partition function as 
a function of the different variables (such as tempera- 
ture, volume, number of particles, etc) amounts to com- 
plete knowledge of the system (with a fixed number of 
particles, since it is the canonical partition function) in 
thermal equilibrium. Note that the computation of the 
partition function is usually regarded as a hard problem, 
since it involves the sum of an exponential number of 
terms. 

To have an intuitive picture of a gauge symmetry, one 
can think of every lattice site as the origin of a coordi- 
nate system [13] . Then, the physics of, say, a rigid rotator 
should be independent of whether it is described from the 
space coordinate system, namely, the coordinate system 
at site i, or from the body coordinate system, namely, 
the coordinate system at site j, where are linked 

by an edge. In this analogy, the gauge field that is at- 
tached, by definition, to the edge (i, j) corresponds to a 
"rotation matrix" in going from one coordinate system 
to the other. Then, rotating the local coordinate system 
at i corresponds to applying a gauge transformation at i, 
which affects all gauge fields at edges adjacent to vertex i. 
Similarly fixing the orientation of the coordinate system 
at i with respect to coordinate system at j corresponds 
to fixing the gauge field at edge (i,j). 

Comparison with an SSM. In order to illustrate what 
a global symmetry is, we consider the Hamilton function 



5 



of an Ising model: 



#1 



Ising 



= - J 



<i,j> 



(6) 



where the spins take values s» = ±1 [note that this is 
different from Eq. (pjj], and the sum is over nearest neigh- 
bors. Here it is clear that a global flip of all spins leaves 
the Hamilton function invariant, whereas any local spin 
flip would amount to some energy change. 



model, and show that a polynomial overhead is re- 
quired in all relevant cases; 



(vii) In Sec. Ill G we show that the completeness of the 
4D Z2 LGT also holds approximately for contin- 
uous models, including all U(l) LGTs as well as 
(Abelian) continuous SSMs; 



(viii) Finally, in Sec. Ill H we discuss our main result in 



terms of the space of all theories. 



III. COMPLETENESS OF THE 4D 



LGT 



A. Quantum formulation 



In this section we prove the main result of this paper, 
namely the completeness of the 4D Z2 LGT. This means 
that the partition function of any Abelian discrete clas- 
sical spin model (including Z 9 LGTs as well as discrete 
SSMs) equals the partition function of an enlarged, in- 
homogeneous 4D Z2 LGT. The couplings of the 4D Z2 
LGT are precisely the parameters that have to be tuned 
so that its partition function equals one specific target 
model and not another. 

As mentioned in the introduction, similar results were 
obtained in [71 H] , but they either relied on having imag- 
inary couplings [7J, or where much more restricted [8]. 
The result presented here is both general and relies on 
real couplings. Although it was essentially given in [3], 
here we will present it in greater detail and rigor, as spec- 
ified in the introduction. 

In order to prove the main result, we will proceed as 
follows: 



(i) In Sec. Ill A will present a quantum formulation of 
the partition function of all Z 9 LGTs; 



(ii) In Sec. 
tion ru 



III B we will present the merge and the dele- 



e, which are the two tools that will enable us 
to transform one interaction pattern into another; 



(iii) In Sec. Ill C we will prove that the Hamilton func- 
tion of an interaction pattern made of all A:-cliques 
of n 2-level particles, for k — 0, 1, . . . , n, (the super- 
clique), with appropriate coupling strengths, equals 
a general Hamilton function on these n particles; 



(iv) In Sec. Ill D we will use the merge and deletion rule, 



as well as the gauge fixing of edges, to construct a 
superclique out of a 4D Z2 LGT; 



Here we present a quantum formulation of the partition 
function ([5]). Similar formulations have been presented 
in pj [20], and have proven to be useful, inasmuch as they 
have used quantum mechanical tools to gain insight into 
SSMs. The basic idea is that manipulations of the quan- 
tum states that appear in the quantum formulation of 
a model allow one to transform them into the quantum 
formulation of another model. Via this detour into quan- 
tum mechanics, one actually maps the partition function 
of the former model into the partition function of the 
latter. 

We define an (unnormalizcd) quantum state of \F\ 
(j— level quantum particles: 



IV' 



LGT/ 



E®i(E- 



e Jmod q) f 



(7) 



where the sum is taken over all spins at the edges e which 
are at the boundary of /, df. That is, one quantum 
system is placed at every face of the lattice in order to 
characterize the interaction at that face (see Fig.[T]). We 
also define an (unnormalized) complete product state 



\a) 



(8) 



with 



!«}/= E ^™#^-^)]|(^ Se)mod? ), (9) 

That is, I a) / contains the factor with which face / con- 
tributes to the partition function; we will sometimes 
write this state as \a(Jf))f to emphasize its dependence 
on Jf . The basis states in ^ and (|9| are the eigenstates 



(v) In Sec. Mm we state the main result of this pa- of the quantum phase shift operator Z\j) := e 



i2nj/q 



per, namely the completeness of the 4D Z 2 LGT. 
This follows from the construction of the super- 
clique (Sec. HID), and from the equivalence of the 
superclique to any Hamilton function (Sec. Ill C ) . 



Ii), 

for j = 0,1,..., q- 1. 

Then, the quantum formulation of the partition func- 
tion of the Abelian discrete LGT ^ is obtained by com- 
puting the scalar product betweenfV'LGT) and |a): 



Then, we will discuss the following aspects of our re- 
sult: 



Zlgt = HV'lgt) 



(10) 



(vi) In Sec. IIIF we investigate how the system size of 



the 4D Z2 LGT scales with the size of the target 



The proof of the previous equality is straightforward, as 
one simply needs to write down explicitly the quantity 
on the right hand side (r.h.s.) of (10) and see that it 



G 



coincides with the partition function Z. Notice that \ip) 
contains the information of the interaction pattern of the 
model (i.e. the lattice on which it is defined), whereas 
| a) encodes the coupling strengths and the temperature 
of the model. Note also that the normalization of the 
states IV'lgt) & n d \a) would amount to an additional 
prefactor in the equality (10), which would not change 



the physics described by that partition function. 
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Figure 1: The state IV'lgt) places one quantum particle (blue 
dots, labeled with numbers) at each face, thereby characteriz- 
ing the interaction of classical spins (black dots, labeled with 
letters) around that face. The quantum spin at face fi, s/ x is 
obtained by summing modulo q the values of the spins at the 
boundary of the face. 



Stabilizer state. 
lizer state EH E21 



The state IV'lgt) [Eq. Q] is a stabi- 
, since it can be rewritten as 



\1> 



LGT, 



(ii) 



where A is a \F\ x \E\ matrix whose entry (/, e) is one if e 
is a boundary edge to the face /, and it is zero otherwise. 
In a d-dimensional square lattice with periodic boundary 
conditions, each face has four boundary edges, and each 
edge is boundary to 2(d— 1) faces. Thus, in this case, each 
row (face) contains four ones and the rest are zeros, and 
each column (edge) contains 2(d — 1) ones. In the case 
of open boundary conditions, the edges at the boundary 
of the lattice are boundary to only (d— 1) faces, whereas 
all faces still have four boundary edges. 

In order to find the generators of the stabilizer of a 
state such as ( 11 ), one does Gaussian elimination to find 



a maximum set of linearly independent columns of A, 
say k. These are used to construct the generators with 
o x operators of the stabilizer. Then one generates n — k 
linearly independent vectors to all the rest (where n is the 
number of particles of the state), which are used to con- 
struct the generators with a z operators of the stabilizer 
group. 



Note that in the case of the state (11), the rank of 
the matrix A depends on the dimension of the square 
lattice d, i.e. this state will be different depending on 
the dimension of the lattice on which it is defined. More 



precisely, one can see that the state IV'lgt) defined on 
a 2D square lattice with open boundary conditions cor- 
responds to the product state |+)®' F ', where |+) is the 
eigenstate of the Pauli matrix a x , with eigenvalue +1, 
(J a:|+) = |+)- On the other hand, IV'lgt) defined on a 
2D square lattice with periodic boundary conditions cor- 
responds to the state \GHZ) = |0)®l F l + |1)®I F I. For 
higher dimensional lattices, e.g. 3D lattices, the state 
IV'lgt) is less trivial. 

Below, we will use the state |^>lgt) defined on some 
lattice as a resource for measurement -based quantum 
computation [53] in order to prove the completeness of 
that LGT on that lattice. The fact that |?/>lgt) de- 
fined on a 2D lattice contains either no entanglement at 
all (open boundary conditions) or a very small amount 
of it (periodic boundary conditions), and thus are use- 
less states from the point of view of measurement -based 
quantum computation, is in agreement with the fact that 
2D Z 2 LGTs are trivial [5] and cannot be complete. 

As noticed in [20], the fact that IV'lgt) is a stabilizer 
state reveals some symmetries in the partition function. 
That is, because \ip) is left invariant under any opera- 
tor s £ S, s\i/}) — \4>), this translates into the following 
invariance in the partition function 



Z = (a\iP) = (a\s\ij). 



(12) 



This implies that there is another set of couplings, de- 
termined by (a'\ = (a\s that yields the same partition 
function as the original set (a\. 



B. Merge and deletion rules 

We now present two rules, the merge and deletion rule, 
which allow us to manipulate the partition function of a 
model and relate it to the partition function of another 
model. The intuitive picture is that the merge rule ap- 
plied to a model with, say, 4 faces, transforms it to a 
model with 3 faces, one being larger, and containing a 6- 
body instead of a 4-body interaction (see Fig.[2|a)). And 
applying the deletion rule to a face amounts to mapping it 
to a model where there is no such face (see Fig.gb)). Al- 
though these rules can be generally defined for Z q LGTs, 
we will henceforth focus on the case Z2, since this is what 
we require for the proof. 

Merge rule. The rule works by setting the coupling 
strength of a face, say <// , to infinity. In order to see its 
effect, we consider ([9]), and we divide each coefficient by 
a factor of e@ Jf , 

\«)f= E e^M^Sese/ 'j-^K J2 s e ) mod a )/. 

si,...,sjj eGdf 

(13) 

Since this is a rescaling of the energy, this does not modify 
the relevant physics that one can derive from the parti- 



tion function. In Eq. (13) it is clear that when Jf — > 00, 
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Figure 2: (a) The merge and (b) deletion rules are tools that 
allow us to map one interaction pattern to another. Blue lines 
indicate faces where there are interactions 



only the coeffcient with 



cos 



( '5~* s e)mod 2 

q eedf 



= 1, 



(14) 



i.e 



(Seed/ s e)mod 2 = remains non-zero. That is, the 
overlap with \ctf(Jf = oo)) becomes a projection onto the 
\0) f state, and imposes the condition (J2 e edf s e)mod 2 = 
on the remaining terms. Due to this condition, one 
of the spins around / is not free anymore, but equals 
the sum of the other k — 1 spins (since it is mod 2) , say 
s b = s a + s c + s d m Fig- [3] This condition is substituted 
in another face where s& participates, e.g. in Fig. |3j the 



face depending on Sh + s g 



Sb becomes Sh + s q 



s e + s a + s c + sa- Thus, this effectively enlarges the face, 
that is, two 4-body Ising-type interactions have become 
one 6-body Ising-type interaction by means of the merge 
rule. Note that this remaining 6-body interaction has 
a coupling strength given by the face which has been 
enlarged (see Fig. [3]) . 



1 • 1 

II • t 

Jf = oc 


• 1 

> -> • 1 


I r 1 

it • < 

4 


• 

c — 



II t 


ib • < 


• i 

4 


• 

,f 3 1 



Figure 3: Merge rule. Setting Jf = 00 sets the condition 
s a + St + s c + sa — 0, and thus, one of the variables becomes 
dependent, say Sb = s a + s c + Sd- This is substituted in the 
right face, which now depends on Sh+s g +s e +s a +s c + Sd, i.e. a 
6-body Ising-type interaction, with the coupling strength of 
the enlarged face, Jb eg h (now Jadcegh)- 

The concatenation of merge rules and the gauge fixing 
of some particles allow us to achieve fc-body Ising inter- 



actions, for any k. For example, in order to generate a 
5-body interaction, we would apply the same process as 
in Fig. [3j and we would gauge fix one of the spins at the 
boundary, say s a . 

Note that selecting what particle on the boundary of 
/ is dependent on the others is an arbitrary choice. That 
is, the face / in Fig. [3] could have been merged with the 
lower or right face (s c , Sb dependent, respectively; see 
Fig. [4]), or with other faces if had more than 2 dimensions. 
However, all choices yield equivalent partition functions. 
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Figure 4: Once we set Jf = 00 it is arbitrary to choose in 
what direction this face is merged. Here, s c is chosen to be 
the dependent variable, and thus / is merged downwards. 

We also remark that using the merge rule to transform 
a fc-body interaction to a fc'-body interaction, with k' > 
k, is only possible if k > 3, since 



k' = 2k - 2, 



(15) 



(In the case of k = 2 the spins would we sitting in 
the vertices and the interactions would be through the 
edges; however, the argument still holds true: applying 
the merge rule along an edge simply creates more 2-body 
interactions) . 

Deletion rule. This rule is obtained by setting Jf = 0, 
that is, by deleting the interaction at face / (Fig. [5j). 
Note that this corresponds to projecting the face / onto 
the state |+) = |0) + |1), i.e. |a/(J/ = 0)) = |+). 





Figure 5: Deletion rule. Setting J/ = deletes the interaction 
in that face. 



C. Method to obtain general n— body interactions 

Here we show that a totally general interaction be- 
tween n 2-level particles can be generated if all k- 
body Ising-type interactions between these n particles 



are available, for any subset of k particles, and all 
k = 0,1, ... ,n. An interaction pattern between n par- 
ticles with all possible fc-body interactions is called a 
/c-clique. We coin the term superclique for an interaction 
pattern between n particles containing all fc-cliques, for 
all k = 0, 1, 2, . . . , n 24J. A "superclique of Ising-type in- 
teractions" is a superclique such that all its interactions 
are Ising-type; in this work, when we refer to a super- 
clique, we will mean this kind of superclique. Hence, 
we claim that a totally general interaction between n 
particles can be generated by preparing a superclique of 
Ising-type interactions among them, and by tuning its 
coupling strengths appropriately. 

In order to prove the claim, first note that a general 
interaction between n spins corresponds to assigning a 
different energy A s to each spin configuration s. Let us 
indicate with a subindex on the coupling strength which 
particles participate in a given interaction of the super- 
clique; e.g. J123 is the coupling strength of the 3-body 
interaction between si,S2 and S3. Hence, we need to 
show that the coupling strengths in the superclique can 
always be tuned so that 

M-i)° + 

Ji(-i) si + ... + j„(-i) s - + 
Ji 2 (-i) sl+S2 + • • • + + 

Jl2 3 (-l) Sl+S2 + S3 + ■ • • + Jn-2,n-lA-lY n - 2+Sn - 1+Sn ' 



hi.. .„(-!) 



si+...+s„ 



= A S 



is satisfied for arbitrary A s and for all s. 

We remark that the number of parameters is commen- 
surate since there are as many A's as spin configurations, 
thus 2 n , and as many J's as interactions in the super- 
clique, 



k=0 



2 n 



(17) 



Note that this includes a "zero-body interaction" , which 
is a global factor Jo that corresponds to a shift of all 
energies, and thus does not change the physics of the 
model. Since such a "zero-body interaction" cannot be 
prepared in the superclique, we will obtain the partition 
function of the final model up to this factor. 

By defining J as a column vector with all coupling 
strengths, 

J = (^0i hi • • - > 4i hit ■ • • j Jn-l,n: ^123j • • • 1 h...n) T \ 

(18) 

and A as a column vector with one energy A s for each 
column, 



A — (A(o,...,q), A(o,...,i), • • • , X(i,...,i)) rj 



we can rewrite condition ( 16 ) as 



CJ = A, 



(19) 



(20) 



where C is a square matrix with the coefficients of equa- 
tion (16). That is, every row of C is made of the 
factors]- l) s <>, (-l) Sl , . . . , (-l) s ", • • • , (-l) s i+-+ s ™, and 
there is one row for each spin configuration (e.g. the first 
row corresponds to s — (0, 0, . . . , 0) and thus contains 
(-1)°, (-1)°, . . . , (-1)°, • ■ • , (-1)°; the second row cor- 
responds to s = (0, 0, . . . , 0, 1), and so on). 

Thus, we need to show that C can be inverted, i.e. that 
one can always find J as a function of the given energies 
A, 



J = c- x x. 



(21) 



In the following we will show that C cannot only be in- 
verted, but it is also an orthogonal matrix, i.e. C~ l cx 
C T . 

Let us denote by r m (r n ) the mth (nth) row of C. 
Then we want to prove that 



r„ = 



Vn, 



(22) 



where • denotes scalar product. Let r % m be the ith entry 
of that row (which corresponds to a given interaction 1). 
We define S as the set of interactions i for which rows 
r m and r„ have the same sign, 



S := : rl n = 



(23) 



and D as the set of interactions i for which the two rows 
have opposite sign: 



(16) 



D ■= {* : rl ± <}. 



(24) 



It is easy to see that Eq. ( 22 ) is only satisfied if 

\S\ = \D\, (25) 

where 1 1 denotes cardinality. In order to prove the latter 
statement, we will show that each element in S can be 
paired with one element in D. 

Now, recall that row r m is obtained by substituting 
the value of the spin configuration, say, s m into the co- 
efficients of Eq. (16), and similarly for row r„ with the 
spin configuration, say, s n . If we compare the spin con- 
figurations s m and s„ they will differ in some positions 
(at least 1 and at most n). Let us denote by X the set of 
spins which have a different value in s„ when compared 
with s m . We pick one element of X which we denote 
by x (that is, x is one spin, for example S3 which, e.g., 
in s„ it takes the value and in s m it takes the value 
1). If a given interaction i contains x (e.g., the interac- 
tion si + S2 + S3 contains S3), then we pair it with an 
interaction i' which equals i but does not contain x (e.g. 
with si + S2). And, conversely, if i does not contain x 
(e.g. Si + s 2 ), then we pair it with i' which equals i but 
also contains x (e.g. s± + s 2 + S3). 

Let us suppose that interaction i belongs to S. Then, 
interaction i' has either added or removed one element 
of X, that is, one element that has a different sign in s„ 
compared to s m . Thus i' will belong to D. Similarly, if 



i initially belonged to D, i' will belong to S. Since this 
argument holds for all interactions i, we have proved that 
we can pair each element in S with one unique element in 
D. Thus, the two sets must have equal size, viz. Eq. (25) 
holds. 

Since the above argument holds for any pair of rows 
r„,r TO , this shows that Eq. (22 1 holds and that C is an 



orthogonal matrix. This concludes the proof. 

Let us remark here that in a superclique each particle 
participates in all possible fc-body interactions, for any 
k = 0, 1, . . . , n, that is, in 



E 

k=0 



n — 1 
k 



2 n- 



(26) 



interactions [25] . This fact will be important in Sec. Ill D 
where it will be the reason to require a 4D instead of a 
3D lattice. 

We stress that this mapping is at the level of the Hamil- 
ton function, in contrast with the other mappings used 
for the completeness results, which are at the level of the 
partition function. In plain words, we have proven that 
the Hamilton function of a totally general interaction 
of n 2-level particles (e.g. including complicated many- 
body interactions, etc) equals the Hamilton function of 
a complicated interaction pattern (a superclique), but 
with simple interactions (Ising-type interactions). This 
mapping may have applications in the Hamiltonian for- 
mulation of LGTs |13) and in their renormalization group 
analysis [25] . 



D. Explicit construction of the superclique 

In the following we show that we can construct a su- 
perclique of Ising-type interactions from a 4D Z 2 LGT. 
Because of the result of the previous section, this means 
that, by tuning the coupling strengths of the Ising-type 
interactions in the superclique, this model can specialize 
to any other (Abelian, discrete) classical spin model. 

In order to generate the superclique starting from the 
4D Z 2 LGT we will only make use of the merge and 
deletion rule, and of the gauge fixing. 

We will first show the generation of fc-body Ising-type 
interactions, for any k = 1, . . . , n in a 3D Z 2 LGT. Then 
we will argue that the fourth dimension is needed to repli- 
cate the spins, and thereby let each spin participate in 
all the interactions required in the superclique. 

First, a "single-body" Ising-type interaction of S\ 
(analogous to a magnetic field), i.e. J\{— is obtained 
by letting s\ interact with all other spins around a face 
fixed by the gauge (see Fig. [6]ja)). 

A 2-body Ising-type interaction is obtained by con- 
catenating the merge rule on the front, lower and back 
face of a cube and creating the face with blue bound- 
aries of Fig. [6jb) . This face depends on all of the spins 
at its boundary (i.e. all spins which are attached to the 
thick, blue line in Fig.[6]jb), e.g. the upper, right, and left 



ji 



J~1 s 3 



Figure 6: "Logical" spins (i.e. the spins that will participate in 
the final superclique) are marked in bold black. Blue shaded 
faces indicate merged faces as in Fig. [3] and spins fixed by 
the gauge are marked in red. Figures (a), (b) and (c) show 
a single-body, a 2-body and a 3-body Ising-type interaction 
with coupling strengths Ji, J12 and J123 respectively. Spin r 
does not participate in the interaction because the blue face 
depends on it twice, i.e. r + r — 0, and the same holds for n 
and ra. 



spins on the front face, the lateral spins on the lower face, 
etc). However, six of these spins are fixed by the gauge 
(red spins in Fig. [6jb) ) , that is, their value is fixed to 
zero. Hence, the big blue face only depends on the spins 
which are not fixed, that is on Si + s 2 + r + r = s% + s 2 
(since the sum is performed mod 2). This corresponds 
to the 2-body Ising-type interaction between si and s 2 
with coupling strength determined by the only face which 
has not been merged (the upper face), viz. Ji 2 (— l) Sl+S2 . 
Furthermore, notice that by setting Ji 2 = 00 as well, one 
enforces si + s 2 = 0, i.e. Si — s 2 . This can be seen as a 
"propagation" of the value of si into s 2 . A concatenated 
application of this 2-body interaction results in an effec- 
tive propagation of a spin through a certain path in the 
lattice (see Fig. [7]) . The direction of this propagation can 
be changed by merging all "covering" faces between the 
incoming and the outgoing spin, as indicated in Fig. [8] 
This will be important in the construction of the super- 
clique, where one needs to propagate logical particles in 
the 4D lattice to bring to the place where the interaction 
occurs. 



ril 


f > 











Figure 7: Propagation of the spin si to the spin S3 by concate- 
nating 2-body interactions with J12 = 00. Particle s\ has two 
ends (i.e. faces that can participate in a fc-body interaction): 
itself, and the right face where S3 is. 



In order to generate a 3-body Ising-type interaction, 
one propagates three "logical" spins in the lattice in order 
to bring as close to each other as possible, with the con- 
dition that their "red u-shapes" are not adjacent. Then 
ones merges all but one of the "cover" faces into a large 
blue face as indicated in Fig. IBTc). This blue face now 
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t 




(a) 



Figure 8: (a) Turn in the path from a left -to-right prop- 
agation ("incoming" spin si) to a down-to-up propagation 
("outgoing" spin S2). (b) Similar turn, but here S2 propa- 
gates from back to front. The dependence of each blue face 
on r cancels because it depends twice on it. 



contains the interaction Ji23(-T) Sl+S2+S3 , that is, a 3- 
body Ising-type interaction between Si,S2 and S3 as re- 
quired. The interaction strength J123 is determined by 
the only face that has not been merged, and it is this 
parameter that will be tuned so that the superclique spe- 
cializes to a general n-body interaction. Note that here 
the dependence on r\ and also cancels, since the large 
blue face depends on each of them twice (and the sum is 
mod 2, thus r\ + r\ = 0, and similarly for r^. 

4-body interactions can be generated similarly. In 
this case, the four logical particles are distributed close 
to each other as shown in Fig. [9j and all cover faces 
are merged to give rise to the Ising-type interaction 
^i234(-l) Sl+S2+S3+S4 - The procedure for the 5-body in- 
teraction analogous; in this case, one adds one more log- 
ical spin and merges the overall cover face (see 10). 




Figure 9: 4-body Ising-type interaction between spins 
si,S2,S3 and S4 with coupling strength J1234. See the cap- 
tion of Fig. [6] for an explanation of the symbols. 

The generalization to fc-body Ising-type interactions, 
for arbitrary fc, is straightforward. First, one propagates 
each of the "logical spins" si,...,s/- in the lattice un- 
til they are as close to each other as possible, forming 
a "rectangular" shape without their red u-shapes touch- 
ing each other (as is the case for the 2-, 3-, 4- and 5- 
body interaction shown above). One can imagine this as 
adding more logical particles on the right of the 5-body 
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Figure 10: 5-body Ising-type interaction 

Ji234 5 (-I) si+S2+S3+Si+S5 ■ See the caption of Fig. § for 
an explanation of the symbols. 



in length, until one has fc logical particles, analogously 
to how particles have been "added" in the generation 
of the 5-body interaction when compared with a 3- or 
4-body interaction. Then one merges all cover faces ex- 
cept for one. This renders the fc-body Ising-type interac- 
tion J\...k{— l) Sl+ ■■■ +Sk , where the coupling J\...k is deter- 
mined by the only cover face that has not been merged. 
The dependence on all auxiliary spins r 1: r 2 , . . . will can- 
cel because, by construction, the boundary of the merged 
face depends on them twice. The coupling strengths J\...k 
will be tuned so that the Hamilton function of the super- 
clique equals the Hamilton function of the specific final 
model (which we will refer to as the "target" model), 
according to Eq. (21). 

Thus, we have shown how to obtain fc— body Ising- 
type interactions, for any k = 1, . . . , n. Now we must let 
each spin participate in 2™" 1 interactions, as pointed out 
in (26). However, we have seen that a spin propagates 



interaction of Fig. 10 thereby enlarging the "rectangle" 



as in Fig. [7J and this propagation ends in a certain face 
(called an "end") that participates in a fc-body interac- 
tion. There it is clear that spin Si can only participate in 
two interactions, corresponding to the left and right ends. 
More generally, the number of ends that an object (or en- 
coded particle) of dimension d e in a lattice of dimension 
d has are 2{d— d e ). Here the logical spin is never propa- 
gated alone, but always "carries" the other three spins of 
an adjacent face fixed (i.e. the shape of Fig.|6ja) is prop- 
agated), hence we essentially have d e — 2. So, for d = 3 
the particle is blocked to have only 2 ends. We need to 
resort to a 4D lattice to obtain 2(d — d e ) > 2 ends (see 
Fig. [IT] for a replication in four dimensions of one spin 
into five other ends). Then, this replication procedure 
can be multiply applied until the particle has 2™ _1 ends, 
that is, one end for each interaction. Note that in this 
replication procedure no loops of spins fixed by the gauge 
are formed. 

We remark that all faces which are not mentioned in 
this construction have to be deleted using the deletion 
rule. We also mention that we have tried several other 
procedures in order to obtain this result in 3D, but none 
of them could avoid the formation of loops of edges fixed 
by the gauge. 

The specific layout of interactions in the superclique is 
the following. The logical particles are distributed along 
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Figure 12: 3D space with w = 0. Logical particles are dis- 
tributed along the x direction and they are propagated along 
the y direction. 



Figure 11: Replication of spins in a 4D lattice: si is replicated 
into S3, S5, S6, ss, sg. Yellow faces have the same meaning as 
blue faces, that is, S2 propagates into S3 by the same means 
as it propagates into S7. Note that no loops of red spins are 
formed. 



the x direction with one "idle" space among them (see 
Fig. 15). Then each of them is propagated in the y di- 
rection. The idea is to use this 3D space (i.e. the space 
with w = 0) to propagate the particles, and to use the 
3D space defined by w = 1 to create the interactions 



required for the supercliquc. 



For example, in Fig. 12 we see 3 logical particles along 
the x direction which are propagated in the y direction. 
Then, at some sites, they are also propagated to the space 
defined by w = 1. In particular, in the first site, they are 
all propagated to w — 1, where the 1-body interaction of 
each of the particles will take place. In the following site, 
Si and S2 are propagated to w = 1, in order to gener- 
ate the 2-body interaction among them. After some idle 
propagation in the y direction (precisely, A(2) idle sites, 
as will be explained below), s\ and S3 are propagated 
to the space w = 1, where they will generate a 2-body 
interaction among them. This goes on for all 2-body in- 
teractions, then for all 3-body interactions, 4-body, and 
so on, up to the n-body interaction. For example, the 4- 
body interaction between si, S2, S3, S4 and then between 



s i, s 2, S3, S5 are shown in Figs. 13 14 



Note that the idle space one has to leave in the y direc- 
tion among each propagation to the w = 1 space depends 
on k. This is because, when the k particles are propa- 
gated to the w = 1 space, they are distributed in a line. 
There, one has to rearrange them in the rectangular form 
explained in Figs. |9j [10] It follows from the construction 
that this rearrangement requires to leave space 




Figure 13: A 4-body interaction in the superclique. The four 
particles si, S2, S3 and S4 are propagated into the space with 
w = 1, shown here. Here they interact in a 4-body Ising- 
type interaction as the one presented in Fig. [9] with coupling 
strength Jijki- Black arrows indicate propagation of the spin 
(as in Figs. [7| [§] the corresponding merged faces are not de- 
picted to avoid overloading). 



between interactions. An overall layout of the propaga- 
tion of particles in the w = space is shown in Fig. 15 



Finally, note that, as indicated in Fig. 15 
a 4D lattice of size 



one requires 



(x,y,Z,w) 



(2n,^A(k) 

k=0 



,1,1) 



(28) 



to generate a superclique of n particles. There is an ex- 
ponential overhead in the system size since one has to 
generate an exponential number of interactions. We re- 
mark that efficient constructions can be found for specific 



A(k) = 2\k/4\ +2 



(27) 



target models, e.g. in Sec. VB wc show that the construc- 
tion of the 2D Ising model only requires a linear overhead. 
We also point out that the construction of the 4-clique 
(i.e. the part of the superclique with 4-body interactions) 
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-4(4) 

Figure 14: View of part of the iv = 1 space. First, parti- 
cles si, S2, S3, S4 are propagated into this space, and they are 
brought close to each other (propagations indicated in black 
arrows) in order to interact in a 4-body interaction with in- 
teraction strength J1234. After A(4) idle particles in the y 
direction, particles si, S2, S3, S5 are propagated into this space 
and, again, they are brought close to each other to interact 
in the 4-body interaction with strength J1235. 











1/ 










Figure 15: 3D space with w = 0. Logical particles are dis- 
tributed along the x direction and they are propagated along 
the y direction. A detailed part of this space is shown in 
Fig. [12] 

is the essential ingredient of the mean-field theory that 
we will construct in Sec. IVIB1 



E. Main result 

Now we can finally gather the results we have proven 
in the last sections in order to state our main result. In 
Sec. |IIID| we have shown that by setting some coupling 
strengths to infinity or zero (merge or deletion rule), the 
partition function of a 4D I2 LGT can become the par- 
tition function of a superclique. Then, in Sec. |III C| we 
have shown that if one tunes the coupling strengths of 
a superclique appropriately, its Hamilton function spe- 
cializes to a totally general Hamilton function between 
n 2-level particles, and thus the corresponding partition 
functions are also equal. Therefore, we have shown that 
the partition function of an enlarged 4D Z2 LGT with 
appropriate inhomogeneous coupling strengths can spe- 
cialize to the partition function of any Hamilton function 
of n 2-level particles. More specifically, we have seen 



that, for any classical spin system, there is a subsystem 
of the complete model (the superclique) that behaves like 
it (when the appropriate coupling strengths are set on it). 

Let us elaborate on the class of models that are em- 
braced by this result. First of all, the completeness result 
holds for models with an arbitrary interaction pattern be- 
tween these n 2-level particles, which includes 

• Models in regular lattices in arbitrary dimension; 
e.g. an 8D Z2 LGT can be mapped to an enlarged 
4D Z2 LGT with appropriate inhomogeneous cou- 
plings; 

• Models on arbitrary graphs; e.g. a Z2 LGT defined 
on a complicated, irregular graph (note that usually 
Abelian discrete LGTs are only defined on hyper- 
cubic lattices and here a much more general class 
is considered); 

• Models with different number of particles partici- 
pating in the many-body interactions; e.g. models 
with 6-body interactions can be mapped to the 4D 
Z2 LGT, which has 4-body interactions. 

• Models with different types of many-body interac- 
tions; e.g. models containing more general 4-body 
interactions (the most general case being to assign 
a different energy to each of the 16 configurations of 
the 4 particles) can be mapped to the 4D Z2 LGT, 
which only contains Ising-type interactions. 

In the second place, notice that the information of 
whether the model possesses a global or a local symme- 
try is also encoded in the Hamilton function. Since all 
Hamilton functions are included in our result, this means 
that the completeness result is valid for 

• Models with local symmetries, i.e. other Abelian 
discrete LGTs; 

• Models with global symmetries, i.e. Abelian dis- 
crete SSMs; e.g. the Ising model or the Potts model 
can be mapped to a model with local symmetries, 
the 4D Z 2 LGT. 

We will discuss how models with different types of sym- 
metries can be mapped to each other in Sec. |VB| 

Furthermore, the completeness result also includes 
general Hamilton functions between q— level particles, 
since one just needs to encode each g-level particle in 
uiq = [log q] 2-level particles. Then, a totally general in- 
teraction between n' q— level particles is generated with a 
superclique of n 2-level particles, with n = n'm q . Thus, 
our result also holds for 

• Models whose particles have an arbitrary number of 
levels; e.g. Z 9 LGTs, which have g-level particles, 
can be mapped to the 4D Z2 LGT, which has 2- 
level particles. 

In conclusion, we have shown that the 4D Z 2 LGT 
is complete for all Abelian discrete classical spin models, 
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including all Abelian discrete LGTs and Abelian discrete 
SSMs. In symbols, the main result of this work can be 
summarized as 



^Abc 



discrete classical 



(J) 



^4DZ 2 LGT 



(J, J') (29) 



where J is the set of couplings in the target model, and 
J' is the set of couplings in the additional particles of the 
complete model. 



F. Efficiency results 

We have emphasized that all completeness results re- 
quire a larger, inhomogeneous complete model when com- 
pared to the target model. Here we investigate how the 
number of particles in the complete model n! scales with 
the number of particles of the target model n. That is, we 
study how the system size of the complete model increase 
when the system size of the target model increases. 

First, we focus on the number of particles participat- 
ing in interactions in the target Hamilton function. If 
the target Hamilton function contains at most k— body 
interactions (and q = 2), in general one needs to generate 
a superclique of k particles in the 4D Z 2 LGT for each 



of these interactions (because the method of Sec. IIIC 



could be applied in this case). This superclique contains 
2 k Ising-type interactions, and thus the same order (up 
to polynomial factors) of particles in the complete model. 
If the target Hamilton function contains M such terms, 
we need to generate each of these interactions, and hence 
require a scaling poly(M, 2 k ). 

However, we notice that the number of particles par- 
ticipating in an interaction, k, usually does not grow with 
the system size. This is the case for basically all SSMs, 
such as the Ising, Potts or clock model, as well as vertex 
models (such as the 6-vertex or 8-vertex model). When 
this condition holds, the scaling of the complete model 
with k is constant. We also note that for particular tar- 
get models, there can be more efficient ways to create 
those interactions, and one does not need to create a su- 
perclique. 

Next, we consider the case of g-state models. Each q- 
level particle requires first an encoding into m q = [log 2 q~\ 
2-level particles. Thus, a general fc-body interaction be- 
tween g-level particles requires to prepare a superclique 
of km q particles. Now the same considerations as above 
apply, and thus a Hamilton function of M terms with 
at most fc-body interactions between q— level particles re- 
quires an overhead 



poly(M, 2 k ,q) 



(30) 



in the 4D Z 2 LGT. 

In summary, the target model can be prepared effi- 
ciently if k scales not faster than logarithmically, and q 
and M scale polynomially with the system size. 

These criteria determine whether a given continuous 
SSM can be approximated efficiently. For example, an 



Ising model defined on a lattice has constant q = 2 and 
k = 2 for all sizes, and as the size of the lattice increases, 
the number of terms in its Hamilton function M increases 
polynomially. Hence, it can be prepared efficiently from 



a 4D Z 2 LGT. Even more, in Sec. VB we will present 
an optimized construction of the 2D Ising model which 
scales linearly with the system size. 

Regarding Abelian discrete LGTs, they are usually de- 
fined on a hypercubic lattice, i.e. with k — 4 fixed, and 
with for some particular, fixed value of q. Thus, as we let 
the system size increase, only M increases polynomially, 
and therefore they can be efficiently prepared from the 
4D Z 2 LGT. 



G. Approximate completeness for continuous 
models 



The completeness of the 4D Z 2 LGT can be extended 
in an approximate way to continuous models. That is, 
the partition function of an (Abelian) continuous model 
can be expressed, up to a certain accuracy, as a specific 
instance of the partition function of the 4D Z 2 LGT. This 
can be trivially seen by considering a target model with 
g-level particles, and then letting q —> oo. 

As we have seen in Sec. |IIIF| the overhead in system 
size of the 4D Z 2 LGT will be polynomial as long as q 
increases polynomially with the system size. 

This extension includes continuous models with global 
symmetries (i.e. continuous SSMs). For example, it in- 
cludes the classical Heisenberg model 



Hciscnbcrg — 



(31) 



<i-,3> 



where the sum is over nearest neighbors and the spins are 
vectors in a unit 3-dimcnsional sphere, Si £ M 3 , \si\ = 1. 

The result also includes continuous models with local 
symmetries, i.e. U(l) LGTs, which are defined by the 
Hamilton function 



Hu(l) LGT - Jf 

i 



cos 



e&df 



(32) 



where the sum is over all faces /, and the cosine of the 
sum of the spins around each face is taken. Here the spin 
variables are angles 9 e £ [0, 2ir). The exploration of exact 
completeness results with continuous models is ongoing 
work. 



H. Space of all theories 

In the previous sections we have established a relation 
that embraces models with very different features, such 
as different dimensions, number of levels of each particle 
or different types of symmetries. The purpose of this 
section is to picture this result in the space of all theories 
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(to be denned below) and to compare it with other similar 
relations that had been obtained previously. 

The space of all theories is a space of models where 
the dimension of the model d is depicted in one axis, the 
number of levels of the particles q in the second axis, and 
the number of particles participating in an interaction 
k in the third [26 . Hence, we identify a specific model 
with a point in this space, which is defined by the triplet 
(d, q, k). Notice, however, that a point actually corre- 
sponds to infinitely many models, since there are many 
features (such as the choice of parameters) that are not 
specified in this space, as we will point out below. The 
space of all theories is useful to visualize different univer- 
sality classes, as we will elaborate on below. 

In the space of theories, our result can be illustrated as 
follows: one can map models with any (d, q, k) (with inho- 
mogeneous coupling strengths and arbitrary patterns of 
interactions) to a model with (4, 2, 4) (see Fig. 16 ). More- 
over, we know that this complete model has Z2 gauge 
invariance, inhomogeneous coupling strengths, and more 
particles than the original model (since it also has all 
auxiliary particles used in the construction of the super- 
clique) . 

This has been proven in different steps. We have shown 
that a model with fixed (d, q, k) can become a model with 
(d',q, k) for arbitrary d! by preparing a fc-clique, and it 
can become a model with arbitrary (d', q, k') by prepar- 
ing the superclique (Sec. HID). Then, the superclique is 



shown to be equivalent to a general interaction between 
n 2-level particles (Sec. IIIC). The mapping to models 



with general q' is achieved by encoding these particles 
into 2 level particles. Note that the mapping of a model 
with a certain d (or q or k) to a model with larger d 
(or q or k) (with more particles) is trivial, since one can 
just leave the extra dimensions (or the extra levels or 
particles) empty. 



the dimension of the model is specified, assuming that it 
is defined on a regular lattice), (ii) inhomogeneous cou- 
pling strengths (here the coupling strength is fixed to J; 
in principle one could generalize it by adding one axis 
for every coupling strength Ji, J2, . . .), (iii) whether the 
model has global or local symmetries, which is related to 
the type of interactions in the model (this could only be 
recognized in the universality classes of the space), (iv) 
whether or not all particles of the model have q levels 
(although it is a natural requirement to ask for the same 
number of levels in all particles) , and (v) the number of 
particles of the model is not represented explicitly. Nev- 
ertheless, this space encodes some of the most relevant 
information and is useful to describe part of our main 
results in a graphical way. 

In order to put the present result into context, we 
include a figure of [8], which represents previous com- 
pleteness results in the space of all theories (Fig. [T7| . 
Fig. [rTfa) represents the completeness result of [7] , where 
it was proven that the partition function of all (Abelian, 
discrete) models can be mapped to the partition func- 
tion of the 2D Ising model with inhomogeneous magnetic 
fields and with complex couplings. Fig. [l7lb) represents 
the completeness results obtained in [5], where it was 
proven that the partition function of a model in d di- 
mensions, with at most fc-body interactions and with 
q-level particles can be mapped to the partition function 
of a model in three dimensions, with fc-body interactions 
and g-level particles. Despite being more restricted, the 
latter completeness result had the advantage of requiring 
only real parameters in the complete model. For exam- 
ple, the result of [8] includes the fact that the partition 
function of a d-dimensional Ising model can be mapped 
to the partition function of a 3D Ising model with real 
couplings. 



Superclique 



Encoding 





5-— q 




Figure 16: The main result of this paper seen in the space of 
all theories: models with arbitrary dimension d, g-level parti- 
cles and fc-body interactions can be mapped to the partition 
function of the 4D Z2 LGT with real parameters. The latter 
can specialize to models with larger d or k by constructing the 
superclique, and to models with larger q by encoding g-level 
particles into 2-level particles. 

Note that, since this is only a 3-dimensional space, 
many properties of the models are not captured in this 
space, such as (i) the specific interaction pattern (i.e. only 



Figure 17: Previous completeness results, (a) All models can 
be mapped to the 2D Ising model with inhomogeneous mag- 
netic fields with complex parameters [7]. (b) Models with 
large (d, q, k) can be mapped to a model with (3, q, k) with 
real parameters |5j. 



It is worth to point out the relation between the com- 
pleteness result for SSMs and the existence of a criti- 
cal dimension defined by the mean-field theory in SSMs. 
Namely, it is known that the critical dimension for classi- 
cal spin systems in SSMs is d = 4. For higher dimensions, 
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a mean-field theory captures the universality properties 
of models. Interestingly enough, the dimension obtained 
with completeness results yields d = 3, which is a theory 
not described by mean-field theory. The result of Fig. [16] 
is more complete in the sense that we have reduced or 
transformed all classical spin models into one point of 
the reduced space, as in Fig. 17 'a), but for all theories 
and with real parameters. 

The notion of criticality in LGTs has been less studied 
numerically. It is known that dD 1 q LGTs in a cubic 
lattice, with d < 4, exhibit two phases: confined and 
deconfincd. In the case d > 4, and q > 5, there is an 
additional phase, the Coulomb phase, that lies between 
the other two [27Tl30| . The situation for non-Abelian 
LGTs is far richer but less understood. 

Universality classes. Note that this space of all theo- 
ries contains models belonging to many different univer- 
sality classes. That is, different families of models (each 
defined by a certain range of d, q, k), where every fam- 
ily behaves similarly around the critical point. Each of 
these models is mapped to the complete model with a 
specific set of parameters, and thus each of these families 
is mapped to the complete model with a specific regime 
of parameters. Thus, the phase diagram of the complete 
model (with one coupling strength in each axis) should 
contain the different universality classes, that is, it should 
mimic the space of all theories. 

In particular, this includes models with global symme- 
tries. We discuss in Sec. |VB| how it is possible that a 
model with local symmetries specializes to models with 
global symmetries is discussed in Sec. |VB| 

Note that the result also includes Z q LGTs with "mat- 
ter fields" . These theories are like the pure gauge theories 
that we have considered so far, and they additionally con- 
tain g-level spins at the vertices Sj = 0, 1, . . . , q — 1, for 
every i £ V. These variables transform under the gauge 
transformation as 



gi(Si) = (Si + l)mod 



(33) 



The gauge-invariant Hamilton function of this theory is 
of the form 



H(s) = Jf cos ( ( s ij + s jk + Ski + S/i)mod ■ 



~ Jv J C ° S ( ~(~ S i + S « + S j")mod 9 ] (34) 

where i j is an oriented edge, and i (j) is the tail 
(head), and i,j,k,l are the vertices at the boundary of 
face / (note that we have changed the notation with re- 
spect to Eq. ([3]) for convenience). We remark that this 
type of matter fields do not correspond to fcrmions since, 
in that case, we should introduce Grassmann variables at 
the site, instead of spin variables. In conclusion, we see 
how the partition function of an (enlarged) pure LGT 
(the 4D Z 2 LGT) can also equal the partition function of 
an LGT with matter fields. Physically, this means that 



these matter spins can be "reabsorbed" into a pure LGT 
which is larger in size and with the proper couplings. 



IV. COMPUTING OBSERVABLES FROM THE 
COMPLETE MODEL 

In this section, we want to point out that our result 
should be used with care for the following reason. The 
computation of physical quantities of the target model 
usually requires one to take derivatives of the partition 
function with respect to some variable. When the parti- 
tion function of the target model is expressed as a special 
instance of the partition function of the complete model, 
the complete model contains more particles, and hence 
more variables (e.g. the coupling strengths of these addi- 
tional particles). However, the derivatives in the target 
model must be taken only with respect to the variables 
in the target model. This is because the couplings in 
the additional particles in the complete model have to 
be fixed, since they are used to transform the interac- 
tion pattern of the complete model to that of the target 
model. Therefore, they cannot be regarded as variables 
from the physical point of view in the target model. 

For example, the mean energy of a (target) model is 
obtained as [3] 



Ul.nr — 



d_ 

"dp 



lnZ tar (J) 



(35) 



We can rewrite this expression making use of our map- 
ping, 



Z tar (/3J) = Z com ( J, J') = J2 



■(J)~P'H add (.J') 



where J is the set of couplings in the target model, and 
J' are the couplings for the additional particles of the 
complete model. Then, we can use (36) to rewrite (35) 

as 



U iax = —^laZ com (0J,p'j'), 



(37) 



where the derivative has to be taken only with respect to 
j3, and not with respect to j3' . 

Other examples of thermodynamical quantities than 
can be obtained from the partition function are the 
Helmholtz free energy [3] 



and derivatives thereof, such as the entropy S, 

'OA" 



S 



or the chemical potential fi 



dT 



V,N 



OA 
ON 



(38) 



(39) 



(40) 



T,V 
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In all of these cases the derivatives should be taken only 
with respect to the variables in the target model, not 
including the additional particles of the complete model. 

In other words, there is a subsystem of the complete 
model which behaves as a physical system, in the sense 
that one can modify, say, its temperature and its physical 
properties vary. More precisely, this subsystem behaves 
like the target model. However, the rest of the complete 
model has a different, auxiliary status: it does not behave 
as a usual physical system, but its couplings and temper- 
ature must have fixed values, which are chosen such that 
the subsystem behaves like the target model. 



V. APPLICATIONS: EXPLICIT 
CALCULATIONS 

In this section we illustrate some of the uses of our main 
result by showing how to explicitly compute quantities of 
the target model as a function of the partition function 
of the complete model. In particular, in Sec. |V A| we will 
show how to derive the order parameter of LGTs, the 
Wilson loop, from the partition function of the complete 
model, whereas in Sec. |VB| we will show how to compute 
the magnetization, which is the order parameter of some 
SSMs such as the Ising model, from the partition function 
of the complete model. In the latter section we will also 
discuss how a model with local symmetries (the 4D Z 2 
LGT) can specialize to a model with global symmetry 
(such as the Ising model). 



where one takes derivatives with respect to all faces en- 
closed by the loop C. 

Here Z sourccs is our target partition function. That 
is, we would first generate a supcrcliquc from the 4D 
Z 2 LGT, and then tune its coupling strengths so that 
its Hamilton function equals the exponent of Eq. ( 43 1 . 



Consequently, the Wilson loop of the target model can 
be computed as 



(W{C)) 



d 
dh 



f 



lnZ com (J, J') 



(45) 



The Wilson loop is also defined for Z q LGTs, and this 
derivation can also be carried over for that case. 



B. Magnetization of the 2D Ising model 

Here we want to illustrate our mapping when we take 
an SSM as our target model. More precisely, we consider 
the 2D Ising model without fields, and we compute its 
magnetization from the partition function of the 4D Z 2 
LGT. The magnetization is defined as 



M tar = 



d_ 
Oh 



LnZ tar (J, h) 



h-yO 



where 



5^ = 



(J)-(V<s)>i£i«4) 



(46) 



(47) 



A. Wilson loops of lattice gauge theories 

A crucial quantity in an LGT is the Wilson loop, which 
is defined as a product of variables in the edges (gauge 
fields) over a closed loop C [2]. In a Z 2 LGT, the Wilson 
loop takes the form 

W(C) = e iw ^^ s ". (41) 



The expectation value of a Wilson loop 

(W(C)) = I^y^-ec^e-^'O, 



(42) 



determines the phase in which the LGT is: in the con- 
fined phase (W(C)} decays as the area enclosed by the 
contour C, whereas (W(C)) decays with the perimeter of 
C in the deconfined phase. This quantity can be derived 
from a partition function with sources [1]: 



J sources 



y^ e -/Ki?(s)-i£ /S F>vs/) 

s 

where sj = e l7r ^c 5 «. Specifically, 
(W(Q) 



(43) 



fee 



d 
dh f 



(44) 



Now, using Eq. ( 29 ) we can rewrite Eq. ( 46 ) as 



d_ 
dh 



lnZ com (J,J',h) 



where 



^com — / c 



t,r(J)+H s . dd (J')-(l/l3)hY, i s i ) 



(48) 



(49) 



Note that here {s^} are the particles in the target model. 
It is known that the magnetization of the 2D Ising model 
is its order parameter, and that it is non-zero, M tar 7^ 0, 
in the so-called ordered phase [3]. On the other hand, 
Elitzur's theorem [TT] states that the magnetization of a 
model with local symmetries must remain zero. 

This applies, in particular, to the magnetization of our 
complete model, the 4D Z 2 LGT. That is, the r.h.s. of 
Eq. (48) corresponds to the magnetization of the com- 
plete model (with couplings J, J'), M com , and this quan- 
tity must be always zero, M com = 0. 

This apparent contradiction is due to the fact that 
Elitzur's Theorem only holds when the couplings are fi- 
nite and, in order to construct the target interaction pat- 
tern starting from the 4D Z 2 LGT, we require to fix cou- 
plings to infinity, i.e. J'j = 00 for some faces /. This 
effectively amounts to breaking the local symmetry on 
that face, since this imposes a particular value on the 
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spins around that face [Eq. (14)]. This explains how the 



partition function of a model with local symmetries can 
specialize to the partition function of a model with global 
symmetries: only if the former model contains "extreme" 
values of the coupling strengths (J 1 * — oo), can the local 
symmetry break spontaneously and behave as a global 
symmetry. 

Explicit construction of the 2D Ising model. In 
Sec. |III D| we have shown the general procedure to con- 
struct an interaction pattern starting from the 4D Z2 
LGT (construction of the superclique) . This procedure 
can be optimized for a specific target model, reducing 
(drastically) the scaling of the size of the 4D Z 2 LGT as 
a function of the size of the target model. In what follows 
we illustrate this by constructing a 2D Ising model from 
the 4D Z2 LGT in a direct way without preparing the 
superclique. 

The Hamilton function of a 2D Ising model is given by 
Eq. ([6]), with the nearest neighbors defining a 2D square 
lattice. Here we will construct this interaction pattern 
directly, without preparing first a superclique. 

First of all, we concentrate on generating a 1- 
dimensional (ID) array of Ising-type interactions. In 
this case, we just need to generate 2-body interactions 
between nearest neighbors distributed on a row. This is 



precisely what we do on Fig. 18 There, spin si^i interacts 



with si ; 2 in a 2-body Ising-type interaction with strength 
Jf 2 , S1.2 interacts with 5^3 with strength 3 , and so on. 
The 2-body interactions are given by Fig. [6](b). 




Figure 18: Generation of a ID array of Ising-type interac- 
tions: Sij interacts with spin Sij+i via a 2-body Ising-type 
interaction with coupling strength Jij ; ; J+ i. Arrows indicate 
propagation of the spin according to Fig. [7] and the cubes 
marked in blue indicate a 2-body interaction according to 
Fig.^b). 



The next step is to create interactions among several 
of these ID arrays. Since we require three dimensions 
to construct the ID array of Fig. [18] we make use of 
the fourth dimension in order to link them as shown in 
Fig. [19} The yellow cubes have the same meaning as 



the blue cubes in Fig. [18] that is, they correspond to 2- 
body Ising-type interactions. In this manner, spin si t i 
interacts with s 2 ,i with strength Ji, 1-2,1, and so on. This 
completes the construction of the 2D Ising model. 




Figure 19: Construction of the 2D Ising model. Each ID array 
interacts with the next one via the fourth dimension, that is, 
Sij interacts with Si+ij via a yellow face, with interaction 
strength J;j ; i+ij. Every layer for different w corresponds to 
the ID array of interactions of Fig. 
shown to avoid overloading). As in Fig. 



(b lue cubes are not 
yellow cubes 



11 



have the same meaning as blue cubes, i.e. 2-body Ising-type 
interactions. 



As can be observed in Fig. [19] the construction of a 
2D Ising model of size nx m requires a 4D lattice of size 
{x, y, z, w) — (2n, 4, 1, to), i.e. the scaling is linear in the 
system size. 



VI. IMPLICATIONS OF THE MAIN RESULT 

In this section we will draw two implications of the 
main result. First, in Sec. |VI A| we will conclude that 
computing the partition function of 4D Z2 LGT is 
^/=P hard; that is, computationally difficult. Then, in 
Sec. |VIB| we will argue that our result provides a new 
method to compute the mean-field-theory of a Z2 LGT, 
which works for finite dimension. 



Computational complexity of 3D and 4D Z2 
LGT 



Our main result implies, in particular, that the parti- 
tion function of the 2D Ising model with magnetic fields 
can be expressed as a specific instance of the partition 
function of the 4D Z2 LGT. The computation of the 
partition function of the 2D Ising model with magnetic 
fields is a #P-complete problem [3T] -colloquially speak- 
ing, this means that it is computationally difficult [32] . 
Thus, we conclude that computing the partition function 
of the 4D Z 2 LGT in the real parameter regime is #P- 
hard, i.e. at least as hard as the other problem. In other 
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words, we have proven that one can map all models to a 
model which is hard to solve. 

The construction presented above also gives insight 
into the complexity ol the 3D Z2 LGT. More precisely, 
in Sec. |IIID| we saw that a 3D Z2 LGT can prepare 
models with fc-body Ising-type interactions, for any 
k = 1, . . . , n, as long as as every particle participates in 
at most two interactions (this was the limitation of the 
two ends of Fig. [7] that made us move to the 4D lattice) . 
This implies, in particular, that the 3D Z2 LGT must 
be as hard as any vertex model with q = 2 and fc— body 
Ising-type interactions. 

On the other hand, using a method introduced in |33j , 
one can show that approximating the partition function 
of the 3D Z 2 LGT in a certain complex parameter regime 
with polynomial accuracy is as hard as simulating arbi- 
trary quantum computations, i.e. BQP-complete [53] . 



B. Mean-Field Theory 

The mean-field theory of a model is an approximation 
to that model where the interaction of a variable with its 
neighbors is replaced by an interaction of this variable 
with a mean field. In this manner, the theory is reduced 
to a 1-body problem, which is useful to gain insight into 
a theory that is difficult to solve exactly. Thus, there 
are as many ways to construct a mean-field theory of a 
model as ways to average over the influence of neighbor- 
ing variables over a given variable. 

In SSMs, mean-field theories of SSMs are generally 
easy to construct. For example, in the Ising model, the 
mean field is a mean value of the other spins, which es- 
sentially corresponds to the magnetization of the model. 

However, mean-field theories for LGTs are generally 
hard to construct. This is due to Elitzur's theorem [TT] . 
which asserts that the mean value of every variable is 
always zero. This problem was circumvented using a 
saddle-point approximation with the inverse dimension 
1/d as an expansion parameter [35J |35] . The restoration 
of the gauge symmetry in then nontrivial in this expan- 
sion [4j. 

In the following we argue that our proof of the main 
result yields a new method to compute the mean-field 
theory of a d dimensional I2 LGT. This method works for 
fixed d, with d > 4, and does not break gauge invariance. 

The method is based on the construction of the 4- 
clique, that is, the interaction pattern in which every 
particle interacts with all the rest in all possible 4-body 
interactions. The construction of the 4-clique for the 4D 
Z2 LGT was shown in Sec. |IIID[ where we constructed 
a first a 1-clique, then a 2-clique, and so on for all k- 
cliques, with k = 1, . . . , n. Thus, the 4-clique is just one 
part of that construction. Since in a 4-clique every par- 
ticle participates in all possible 4-body interactions, this 
corresponds to the coarse grained construction charac- 
teristic of a mean-field theory. 



We believe that the explicit construction of the 4- 
clique could be useful for computer simulations of the 
mean field theory of an LGT, and for this reason we pre- 
sented it in detail in Sec. IIIIDI We remark that the con- 
struction of a 4-clique of n particles requires a 4D Z2 
LGT of exact size 



(x,y,z,w) = (2n,4 



1,1) 



(50) 



For example, to construct a 4-clique of 8 particles (70 
interactions) a 4D lattice of size (16, 280, 1, 1) is required. 

We remark that the construction of the fc-clique re- 
quires to use the merge rule, and thus to set some 
coupling strengths to infinity (Sec. HID). As discussed 
in Sec. |VB[ Elitzur's Theorem does not apply in this 
regime, that is, gauge invariance is violated. Thus, in 
order to preserve gauge invariance, we can take a large 
but finite value of J, so that gauge invariance is still pre- 
served. This will lead to an approximate construction 
of the 4-clique and thus to an approximate mean-field 
theory. But this does not represent a problem, since the 
mean field theory is an approximate method itself. So 
our method allows to compute an (approximate) mean- 
field theory without breaking gauge invariance. 

This is natural in the sense that it is at the root of 
the difference between SSMs and LGTs: global versus 
local invariance. An example may illustrate this fact. 
Namely, there is a complicated and indirect way to relate 
the 3D Ising model with the 3D Z2 LGT (the duality 
relation mentioned previously) [2J. Without dwelling on 
the details, it needs a transformation that it is non-local 
in the fields, even though each model is local in the fields. 
Therefore, it is natural that the transformation relating 
two types of different types of models, if exact, must be 
singular in the couplings. This is a way of signaling a sort 
of critical point separating two types of different phases 
or models. Our rules allow to establish a more universal 
and direct relation, as we have shown in the previous 
sections. 

Note that our method requires no expansion in the di- 
mension of the lattice 1/d as in [35] |35J , but it works for 
fixed dimension d = 4. It also trivially works for fixed 
dimension d > 4, where one can construct the 4-clique 
by the same method, and simply not use the extra di- 
mensions. Indeed, our method consists on maximally in- 
creasing the connectivity of each of the spins, which can 
be seen as increasing the dimensionality of the lattice 
(since, e.g. in a d-dimensional hypercubic lattice every 
spin participates in 2c? interactions). Note also that we 
can construct a fc-clique for any fc, and thus our method 
does not only apply to the usual case of a hypercubic lat- 
tice, i.e. fc = 4, but it is more general. Finally, notice that 
our construction of the fc-clique requires a very specific 
gauge fixing (compare it with [35], where the mean-field 
theory of Z2 LGTs is discussed with and without gauge 
fixing) . 

We mention that this construction can be generalized 
to Z q LGTs, since one can also define a merge rule and 
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apply a similar construction. 



VII. EXTENSIONS 

In this section we present some extensions of our main 
result. We will argue that there are two other models (the 
3D Z2 LGT and the 3D Ising model with 3-body interac- 
tions) which are also complete if some specific boundary 
conditions are imposed on them. Finally, we will com- 
pare the 3D and the 4D Z2 LGT and discuss whether a 
reduction of our result to the 3D Z2 LGT is possible. 



Completeness of the 3D Z2 LGT with fixed 
boundary conditions 
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Figure 21: (a) There is a loop of spins fixed by the gauge 
(straight, red lines) in the inner part of the lattice, and there 
are spins fixed by the boundary conditions (wavy, red lines) on 
the boundary of the lattice, (b) The loop can be "deformed" 
until one of its edges is the fixed by the boundary condition, 
thereby opening the loop. 



As mentioned in Sec. Ill D the only obstacle in prov- 
ing the completeness of the 3D Z2 LGT was that every 
superclique constructed from this lattice involved loops 



of spins fixed by the gauge (see Fig. 20), which are not 
allowed, i.e. do not leave the physics of the Hamilton 
function invariant. An easy way to overcome this prob- 
lem is to impose fixed boundary conditions in the 3D lat- 
tice. In this case, some spins are fixed to zero due to the 
boundary conditions, not due to the gauge fixing. Then 
every loop that was formed in the superclique (e.g. the 
loop shown in Fig. [20{h)) can be "opened" by replacing 
one spin fixed by the gauge by one fixed by the boundary 
condition (Fig. [Mb)). 





Figure 20: (a) The replication of spins in a 3D Z2 LGT causes 
loops of spins fixed by the gauge (thick, red line), (b) One 
spin in the loop is fixed by the boundary condition (wavy, red 
line), thereby "opening" the loop of spins fixed by the gauge 
(straight, red lines) and allowing replication without loops in 
3D. 

In fact, one only needs to fix these boundary conditions 
at the outer boundary of the 3D lattice. This is because 
every loop in the inner part of the lattice (e.g. Fig.[2l|a)) 
can be "deformed" until one of its edges becomes the 
edge of the outer boundary, which will open the loop (see 
Fig, pjl^b)). This proves that the completeness results of 
the 4D Z 2 LGT hold as well for the 3D Z 2 LGT with 
fixed boundary conditions. 

Although the bulk physics of the 3D Z 2 LGT with 
boundary conditions is the same as that of the 3D Z2 



LGT in the thermodynamic limit, in general boundary 
effects may give rise to new phenomena. 



B. Completeness of a 3D Ising model with 3-body 
interactions and fixed spins 

We remark that the same completeness results could 
have been obtained with a 3D SSM with 3-body Ising- 
type interactions. This is because one can construct a 
superclique out of the this model, since 

(i) this model possesses 3-body interactions, and thus 
the merge rule increases the number of particles 
participating in an interaction [see Eq. (15l], and 



(ii) it is an SSM, thus the fixed spins can form loops, 
hence, an initial 3D lattice suffices. 

Although this model is simpler than the 4D Z 2 LGT (in 
terms of dimensions and many-body interactions) and it 
exhibits the same completeness property, it presents two 
disadvantages. The first one is that, in an SSM, all spins 
fixed inside the lattice must be fixed by the boundary 
conditions, since they do not possess any gauge sym- 
metry. Therefore every fixed spin should be taken into 
account initially in the Hamilton function, which is im- 
practical and difficult to justify physically. The second 
disadvantage is that it is difficult and to define a regu- 
lar 3D lattice only made of 3-body interactions. Since 
the requirement is that the interactions are among more 
than two particles, instead of a lattice with only 3-body 
interactions, we could consider many regular triangular- 
lattices in two dimensions, connected by parallel lines in 
the third dimension. 

Although an Ising model [as defined in Eq. ([6] defined 
on such a lattice, and with many spins fixed in the middle 
of the lattice, would exhibit the same completeness prop- 
erties, its physical interpretation and applications are not 
as clear-cut as those for the 4D Z2 LGT. 
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VIII. CONCLUSIONS AND OUTLOOK 

We have shown that the partition function of any 
Abelian, discrete classical spin model can be expressed 
as the partition function of an enlarged 4-dimensional 
pure lattice gauge theory, with gauge group Z 2 , the 4D 
Z 2 LGT. The values of the coupling strengths of the 4D 
Z 2 LGT determine to what model its partition function 
specializes. The enlargement of the 4D Z 2 LGT is poly- 
nomial in all relevant cases. We have also proven that the 
partition function of Abelian, continuous classical spin 
models can be expressed approximately as the partition 
function of the 4D Z 2 LGT. 

We have illustrated the use of our result by computing 
relevant quantities of the target model as a function of 
the partition function of the 4D Z 2 LGT, and by giving 
a specific example on how the 4D Z 2 LGT specializes 
to the 2D Ising model. We have also discussed how it is 
possible that a model with gauge symmetry can specialize 
to a model with global symmetry. Our result yields a 
new method to compute the mean-field theory for Z 2 
LGTs on hypercubic lattices of dimension d > 4. It also 
allows us to assert that computing the partition function 
of the 4D Z 2 LGT is #P-hard, that is, computationally 
difficult. Finally, we have shown that the same result can 
also be proven with two other models: the 3D Z 2 LGT 
with fixed boundary conditions, and the 3D Ising model 
with 3-body interactions and fixed spins. 

We believe that our results render further insight into 
the structure of Abelian, classical spin models, as they 
provide an explicit mapping between very different mod- 
els (including e.g. models with different types of symme- 
try). In particular, the fact that the phase diagram of 
the complete model must contain all universality classes, 
in principle could be used (i) as a new tool for the study 
of known universality classes, and (ii) as an approach 
to explore possibly unknown universality classes (since it 
embraces a larger class of models than the ones usually 
considered). We also believe that our method to com- 
pute the mean field theory of Z 2 LGTs can be useful for 
computer simulations. 

There are many possible ways in which one could gen- 
eralize and further explore our results. For example, it 



is an open question whether our result can be proven for 
the 3D Z 2 LGT, instead of the 4D, since, in terms of 
the phase diagram, the 3D Z 2 LGT is the simplest non- 
trivial model, for it exhibits a confined and a deconfincd 
phase (this holds for cubic lattices, i.e. k = 4). How- 
ever, our attempts to reduce the result to 3D have failed 
(since loops were always formed), and we do not know if 
this is due to our approach or if this is not possible. One 
reason why this may not be possible is that in a 3D lat- 
tice there are non-trivial knots, whereas in 4D all knots 
are trivial 39 1. Thus, the loops in 3D could be seen as 
knots, and one may only be able to unknot them using 
the fourth dimension. 

Another possible generalization concerns the possibil- 
ity of mapping all known LGTs (including LGTs with 
randomness, and non-Abelian LGTs) to a certain com- 
plete LGT, from which all the rest can be derived. We 
have already mentioned that non-Abelian gauge groups, 
either discrete or continuous, are more involved and de- 
serve a separate study. We have also mentioned that 
LGTs with fermionic matter should be modeled using 
Grassman variables and that these are beyond the scope 
of the mapping presented here. There is a class of random 
classical models, defined both in SSMs and in abelian 
LGTs, that has recently attracted much interest in the 
context of quantum error correction for topological sta- 
bilizer codes [T71 [37J 13H] ■ We envisage the possibility that 
these models with randomness can also be fitted into a 
completeness construction of the type described here. 
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